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USL Program $1400 has been ‘sritten to provide sonar design engineers 
with a technique for selecting filter bands in passive sonar systems ye 
in response to the requirements of reference (a). For a set of rectan- 
gular filters, the program produces plots of output signal-to-noise 
ratio and standard deviation of random bearing error as funotions of 
mid-band frequency. Each set of plots corresponds to a specified oom- 
bination of target, environmental, and array spectral characteristics. 
In addition, the detection performance of the appropriate Eckart filter, 
which is optimum for the particular input parameters utilized, is included 
as a reference standard. Features of the program design, such as the 
format of input data and the control of output information, have been 
selected for ease of use. R 


PROGRAM DESCRIPTION 


Program S1400 represents a first-out response to the requirements 
of reference (a); all input information, which might be stored in a data 
bank or otherwise computed, must now be provided by the user in the form 
of tables. These tables of source level (LS), noise level (IN), pro- 
pagation loss (NW), array direotivity index (NDI), and hydrophone sen- 
sitivity (NHS) are required in dB versus frequency (in kHz). It is 
assumed that each function is then represented by linear segments con- 
necting the given breakpoints in the semi-log domain, a representation 
commonly used by sonar design engineers. The frequency domain to be 
investigated is selected by specifying values of minimum and maximum 
frequency and the number of equally-spaced intermediate points. An 
interpolated value for each of the five input quantities is then stored 
at every frequency, and input signal and noise are computed by: 


SIGNAL(f) LS(f) -Nw(f) + NHS(f) as EN (1) 
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MOISE(f) = LN(f) -MDI(f) + mHS(f) . (8) 


There are situations of interest in which the representation of the 
input noise field by a combination of LN and NDI is neither practical 
nor realistio. In such cases the input sonar self-noise may be read in 
direotly, using a flag to eliminate the use of the MDI table. 


A set of desirable filter bandwidths (in ooteves) is specified in 
the input data. For each bandwidth, caloulations are performed as the 
filter is moved across the allowable frequeney domain. Values of average 
signal (S), everage noise (Nf), and rendom bearing error (SIGTH) are 
associated with the arithmetic mid-bend frequency (F,) for each position- 


ing of the filter. 53 and W are found by using the trapescidal 
integration rule over the linear signal and noise curves and properly 
normali for the band. SIQOTH is specified by equation (6-36) of 
reference (b) as: 


te TY 


235 
SIGTH = 5.3 Fp) So ass, 1+1 =S 3/") (3) 


where 

C = velocity ef sound in feet/second (input), 
D dipole spacing in feet (input), 

T = integrating time in seconds (input), 

NG. oe ener tent "aia nT enw 
F = bandwidth in Hs, eee 


and 


Ky = velocity gain constant, 


The velocity gain constant is also a function of S$ and N , given in 
equation (5-2) of reference (b) as: 


K, = 6.67( SA) ‘ (4) 
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Each value of random bearing error is compared to a computed gste to 

Whether automatic tracking is feasible at the corresponding 
Si. The tracking gate limit is specified by equations (4-39) and 
(7-2) of reference (b) as: 


4 
Gare = sin“1(2:0°C20 ) (5) 
° 


Since SIGTH represents a standard deviation of bearing error, 
is considered acceptable when 3°SIGTH < GATE , and is so indicated on 
the output plot. 


The response of the Eckart filter for each case, given the spectra 
of input signal and noise, is specified in reference (c) by: 


ECK(f) = STGNAL(?f) 


[worsa(t)] * : 


(3/f) pox 22 the ratio of the output signal (s2/n*) to output noise 


(S/N), both averaged over the entire input frequency range. The 
effective signal bandwidth (BWS) is determined by finding the frequency 
limits (in kis), F1S and P28, which are located at the half-power points 
of the output signal curve. This is considered to be a rough first 
approximation, and an improved method for determining F1S and F28 will be 
incorporated in a second version of this program, according to the 
development in reference (d). For passive detection systenus of the type 
proposed, the corresponding recognition differential (in dB) is given 

in reference (e), equation II, as: 


NRD = -19.7 -5.0 10619 (6°Bws) . (7) 


The signal excess (in dB) which would be available using the particular 
Eokart Response function is given by: 


NXS = (Si) x “WARD . (8) 
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In addition, the signal midfrequency (FMS), which divides the output 
signal ourve into two equal power regions, is found. The values of Fis, 
F2s, FMS, RD, NXS, and Ci) are printed on the output plot of the 


Bokart filter response. 
PROGRAM FEATURES 


A complete set of all plots which can be generated by a single set 
of input information is presented in Figures 1-17. ‘There are a totel of 
13 different plot types, as shown in Table I, with plot types 10 and 11 
repeated for each filter bandwidth selected. 


PLT PLOT 
B PLOT DESCRIPTION INDICATORS 
1 Input Source Level (dB) vs. Log Frequency PLOTLS 
8 Input Noise Level (dB) vs. Log Frequency PLOTLN 
3 Input Propagation Loss (dB) vs. Log Frequency PlOnm 
4 Input Direotivity Index (dB) vs. Log Frequency PIOTDI 
5 Input Hydrophone Sensitivity (dB) vs. Log Frequency PLOTHS 
6 Input Signal and Noise (dB) vs. Log Frequency PLOT™ 
1 Input Signel and Noise (linear) vs. Linear LPLOT™ 
Frequency 
8 Input Signal/Moise (dB) vse. Log Frequency PLOTR 
9 Input Signal/Noise (linear) vs. Linear Frequency LPLOTR 


For each bandwidth: 


10 Averege Signal and Average Noise (dB) vs. Log PLOTAN 
Midfrequency 

u Average Signal/Average Noise (dB) vs. Log PLOTAR 
Midfrequency 

For Bokart filter: 

12 Bokart Filter Response (dB) vs. Log Frequency PLOTEF 

13 Eckart Filter 8 1 Response (dB) vs. Lo uveno: PLOTES 


Table I: List of Plots 
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It is clear that a paremetrio study involving many cases could generate 
@ tremendous number of graphs. ‘Therefore, two progrem features have 
been included to minimise the volume of plotted information. ‘The first 
feature is the use of plot indicators in the iaput data to allow the 
user to seiect only those plots which he needs in each case. Since cases 
will differ by a ohange in one or more of the five primary data tables, 
the corresponding plot indicator must be set nonsero as each table is 
fed in. The remaining eight plot indicators will retain the values set 
in the first oase, unless the user purposely reorders them. The second 
progrem feature is the use of subroutine OPILOG to automatically reduce 
the extent of the frequency axis, in each of the semi-log plots, to the 
minimum required to present the information. This subroutine also draws 
: the logarithmically-spaced tic marks end provides appropriate labelling. 


: Subroutine GRAPH is designed to plot the input break-point tables on 
semi-log graphs, taking into account the shift made by subroutine OPXLOG. 
Appropriate y-axis labelling is also provided by this subroutine. 


Subroutine INTERP is designed to operate on a table of break-point 
values, such as those desoribed above for the five input variables, to 
produce dependent values corresponding to an equally-spaced set of 
independent values. The interpolation equation is controlled by two 
indicators, INDX and INDY, corresponding to plot axes x and y, respeo- 
tively. An indicator set to sero means that the corresponding break- 
point variable is plotted linearly; a value of one indicates a log- 
earithmioc plot. For example, a choice of INDX=1 and INDY=0 is 
needed for the semi-log inpyt tables of this program. After extraction 
by subroutine INTERP, the levels in dB are converted to linear units 
prior to integreting over the frequency domsin. 


The input deok for a typical case is presented in Appendix A and the 
printout for this case is given in Appendix B. This printout corresponds 
to the information contained in Figures 1 - 17. The running tise for 
this case, where the frequency interval of 0.5 to 10.0 kis was covered 
in steps of .1 kHz, was approximately 35 seconds on the Univac 1106. 
Progrem S1400 in its ourrent version is available on the first file of 
OUR tape U455. 


INPUT DECK 


The NAMELIST data format has been selected for this program because 
it presents the greatest flexibility and minimum number of restriotions 
for the user, There are no fixed format or order requirements on the 
data cards within each case structure, except for the beginning and end 
ecards. Minor restrictions are as follows: each record card sust not 
start in column 1 and must end with a comma; data items on a single card 
must be separated by commas; and the name of a variable and its essigned 


| ce eaenemaeeaianems 
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value must appear on the same card. The first case must assign values 
to all of the inputs, with the exception of Plot indicators for undesired 
grephs. Thereafter, it is only necessary to provide, in the data cards 
of a case structure, those values which are to be changed from the pre- 
ceeding case inputs. The single exception to this rule is the setting 
of the five input deta plot indicators, which ere serced between each 
case. The usual procedure, then, would be to set the eppropriate data 
Plot indicator nonsero each time a new table is read in. 


The input information for each case structure is contained in the 
following oard set. 


A. First oard - $INPUT in colums 2-7, 
B. Intermediate cards (in any order) - 


1. Case number: CASE 


The progrem automatioally increments the case nusber, starting 
from the value specified in any input case structure. 


8. In olated table desoriptors: MIN. FMAX 


These three values define the frequency bounds (in kits) and the 
number of frequency points in the interpolated input teble. 
All values must be positive, with FMIN < Fal. 


3. Noise flag: NFLAG 


This is a logical variable whioh controls the use of the NDI 
table, as desoribed in “Progrem Desoription” above. MNMFLAG = 
-TRUE. means that MOISE(f) is read in directly; NPLAG = FALSE, 
means that both LN and NDI tables ere required. 


4. Source level data table: FLOTIS, NLS. FLS LS 
ee nn ee 


There are five input data tables, of whioh this is the first. 
PIOTLS 1s the plot indicator which prodwoes the oorrespoading 
output greph (Figure 1) when set non-sero. NLS specifies the 
number of value pairs in the data table to follow. FL3 end LS 
are the tabulated values of frequency (in kis) and source level 
(in dB), which must be punched in the order of their subscripts, 
As an example, a source level table such as: \ 


FLS (tis) LS(dB) 


t 
.5 90.0 

8 96.8 | 
1.0 101,45 


7. 


would be coded as: 
FLS(1)= .5, .6, 1.0, £8(1) = 90.0, 95.8, 101.45, e 

If more than one card is needed to preseat a list, the integer 
im the parenthesis tells the location in the table of the first 
value to the right ef the equal sign. That is, 

SSSSSBees 18(1) = 90.0, 96.8, 

L3(3) = 101.45, See e 

Noise level data table: PLOTLN, MIN, FIN, LN 


Similar to (4), above; sample plet in Figure 8. 


Propagetion loss date table: PLOTS, ea, Fi, we 
Similar to (4), abeve; sample plet in Figure 3. 


Directivity index date table: PLOTDI, MDI, FWDI, MOI 
Similar to (4), above; semple plot in Figure 4. 
sensitivity data table: P Pus 
Similar to (4), above; sample plot in Figure 5. 
Plot indicators for signs] and noise: 
@. PLOTW - plot type 6; sample plot in Figure 6. 
b. LPLOT - plot type 7; sample plot in Figure 7. 
Plot indicators for signel/noise retio: 
@. PLOTR - plot type 8; sample plot in Figure 8. 
b. LPLOTR - plot type 9; sample plot in Figure 9. 
Plot indicators for average values for each bandwidth: 
@. PLOTAN - plot type 10; sample plots in Figures 10, 12 and 14. 
b. PLOTAR - plot type 11; sample plots in Figures 11, 15, and 15. 
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12. Plot indicators for Eckart information: 
@. PLOTEF - plot type 12; sample plot in Figure 16. 
b. PLOTES - plot type 13; sample plot in Figure 17. 
15. Bandwidth specifications: MIN, MAX, DSL: 


‘ 
MIN and MAX define the bandwidth limite (in ocotaves) to be 
analysed by the program, and DEL determines the difference (in 
cotaves) between successive bandwidths. All three values sust 
be positive, with MIN < MAX. 


14. Values for SIGTH end GATE caloulations: Cc, D, TP: 


C is sound velocity (in feet/second), D is dipole spacing (in 
feet), and T is integration time (in seconds). 


15. on: SKIP: 


A value of SKIP = .TRUE. in a case structure will omit the 
execution of that case and al). successive cases until an input 
with SKIP = .FALSB, . 


16. 8 option: STOP: 


Logical variable STOP = .TRUE. will stop the program after the 
execution of the case; if this is not desired, STOP = .FALSE. . 


C. End card - SEND in colwums 2-5. 
FUTURE PLANS 


At present, plans for the second version of Program S1400 inolude 
the following: (1) applying the approach of reference (d) to the 
selection of the Eckart filter frequency limits; (2) investigating methods 
for storing and/or computing the input data tables to simplify perapstrio 
performance studies; and (3) including the effects of non-linear pro- 
cessing in the beamforming simlation. 


L. T. TEIN 
Research Physicist Mathematician 
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